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The long range properties of the compact U(l) lattice gauge theory with 
the multi-level algorithm 
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The 4D compact U(l) lattice gauge tiieory (LGT) in the confinement phase is studied with the multi-level 
algorithm. The static potential, force and flux-tube profile between two static charges are precisely measured 
fi-om correlation functions involving the Polyakov loop. Universality of the coefficient of the 1/r correction to the 
static potential, known as the Liischer term, and the transversal width of the fiux-tube profile are investigated. 



1. Introduction 

If the confining gauge theory is related to an 
effective bosonic string theory, the long range be- 
havior of the potential between static charges sep- 
arated by distance r is expected to have the form: 



Vir)=ar + fi+^ + 0{^). 



(1) 



Here a is the string tension and /i denotes a con- 
stant. The third term is known as the Liischer 
term PP, where the coefficient, 7 = — 7r((i — 2)/24 
with space-time dimension d, is considered to be 
universal. The effective bosonic string theory also 
predicts that the width of the field energy distri- 
bution of the flux tube diverges logarithmically 
as r ^ cxo 2 . Recent Monte-Carlo simulations 
of various LGTs support Eq. ^ with high ac- 
curacy: the confinement phase of ^2 LGT [3], 
SU(2) LGT and SU(3) LGT 

In this report we investigate the long range 
properties of compact U(l) LGT in the confine- 
ment phase. From the Polyakov loop correlation 
function (PLCF), we extract the potential and 
force to see whether this theory also supports the 
presence of the universal correction to the static 
potential. We then measure the profile of the flux 
tube induced by the PLCF, and see the behavior 
of its width as a function of r. However, due to its 
strong coupling nature, the Monte Carlo simula- 
tion of compact U(l) LTG is numerically difficult, 
when its long range properties are of interest. To 
obtain reliable signals, we apply here the multi- 



level algorithm proposed by Liischer and Weisz 
originally for SU(3) LGT [S], which helps to re- 
duce statistical errors exponentially. 

2. Numerical procedures 

We adopt the terminology of Rcf. . To mea- 
sure the PLCF, (P*P(i?)), where R = r/a, we 
take second-level averages of the two-link cor- 
relators T{m;R;i) = U^{m)Ui{m + R i) with 
m = (ms, mt) as 

T^'^\ms,fat]R;i) = [T(m; R]i)T{ra + l; R;i)],{2) 

which is achieved by updating link variables ex- 
cept for the spatial links at to* = 1, 3, . . . , iVf — 1. 
Nt is the timelike extent of the lattice volume. 
i = 1,2,3 denote directions between the two 
charges. We call this procedure the internal up- 
date. We repeat the internal update until reason- 
ably stable values for T^^^ are obtained. Then the 
PLCF at a spatial site rag is constructed as 

P*P{ms;R; i)^Re T'-^\ms, 1; R; i) 
xT(2) (m„ 3; i?; z) • • • T'^) (m„ AT, - 1; R; i). (3) 

The average with respect to rUs and i provides 
[P*P{R)]i^. The desired expectation value is the 
average of [P*P{R)]i^ for = 1,2,..., N^. In 
the actual measurements, we have also applied 
the multi-hit technique to the timelike link vari- 
ables for R> 2 before constructing T. The static 
potential and the corresponding force are taken 
as (neglecting terms of 0{e^'^-^^*)) 
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ay(i?) = -i-ln(F*P(i?)), 



(4) 
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a^F{R) = aV{R) - aV{R - 1) , (5) 

where i? = 1/2 + 0(a). 

In order to measure the flux-tube profile, one 
needs to compute a correlation function of the 
type 

where O is a local operator, (• • ■)j denotes an av- 
erage in the vacuum with the PLCF, and (• • •)o 
an average in the vacuum without such a source. 
To measure {P*PO)o on the mid-plane between 
two static charges, we parameterize the position 
of the local operator n as n — {ni —mi+R/2, rij = 
nij+x, nk = mk+y-, rit =mt) and take the average 
of the two-link-local-operator correlator 

TO^^\ms,mt;n;R;i) = 

[T(m; R; i)0{m + 4; n; R; i)] . (7) 

with 0{m;n;R;i) = U^{m)UiiTn + Ri)0{n) in 
addition to T^^). Combining TO*^^ and T^^), we 
obtain {P*PO)a. 

As local operators, we have used OE{n) = 
'iOf^i^{n) for the electric field, and Ok{n) — 
27rifcp(n) for the monopolc current as defined in 
Ref. |3- Note that the second term of Eq. © 
vanishes for these operators since Oe and Ok are 
parity odd. 
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Figure 1. The static potential as a function of 
r/ro- The constant fi is subtracted by using the fit 
result of Eq. (1), neglecting 0(l/r^) corrections. 
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Figure 2. The force as a function of r/rQ. The 
dotted line corresponds to the theoretical predic- 
tion of Eq. (1) up to 0(l/r^) corrections. 



3. Numerical results 

We generate a sequence of independent gauge 
field configurations, by using the Wilson gauge 
action of compact U(l) LGT on a 16^ lattice at 
f3 = 0.98, 0.99, 1.00, 1.005 and 1.01. For the first 
three /3 values, the configuration has been gen- 
erated after 500 thermalization sweeps, and they 
were separated by 100 sweeps, where one Monte 
Carlo update has been achieved by 1HB/30R. 
For f3 = 1.005 and 1.01 the thermalization sweeps 
have been taken 1000 and 3000, respectively, with 
1HB/50R Monte Carlo update. 

For the static potential and force, we have used 
TVc = 1050, 1250, 2050, 2400 and 3200 configu- 
rations, respectively. The internal updates have 
been taken as iViupd =10000, 8000, 5000, 3000 



and 1000. These numbers are optimized to mea- 
sure the PLCF up to i? = 6 with < 10 % er- 
rors, where the PLCF themselves take values 
from 10^'' to 10^^^. From the force, we have 
introduced a scale based on Sommer's relation 
rQi^(ro) — 1.65. Lattice spacings in units of rg 
are found to be a/ro = 0.470, 0.422, 0.353, 0.305 
and 0.217, respectively. 

In Figs. Hand 121 we show the static potential 
and force from all (3 values. They show good scal- 
ing behaviors. We also plot the theoretical pre- 
diction for the force based on Eq. Q in Fig. |21 
We find that the long distance behavior is well- 
described by the function, F — dV/dr — a — j/r'^ 
with oTg = 1.65 — 7 = 1.39, which contains no 
fitting parameter. This supports the universal- 
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Figure 3. The profiles of the electric field (Ez) 
and monopole current (k,^) for R — 5 a.t (3 = 0.98 
as a function of the radius of the flux tube p. 



ity of 7/r correction to the static potential. Sur- 
prisingly, another feature in common with non- 
Abelian gauge theories is that this function also 
fits the data down to relatively short distance to 
r/rg ~ 0.3. For this, there is as yet no explana- 
tion. 

We have measured flux-tube profiles with 
lengths i? = 3 to 6 at /3 = 0.99 to 1.01 and 
i? = 3 to 5 at /? = 0.98. We have applied 
iViupd = 200 ~ 8000 depending on R and /3. The 
number of configurations are Nc = 300 for all (3 
values. As an example, we show the profiles of the 
electric field and monopole current for i? = 5 at 
/3 = 0.98 in Fig.O which is the longest flux tube 
in these measurements: r/rg = 2.35. Here we 

.y2 



have taken the cylindrical average: p = -^/a 
and ip = tan~^(y/x). As shown in this figure, we 
have obtained very clean signals, however, rota- 
tional invariance for the monopole current profile 
seems not good due to the relatively large lattice 
spacing a.t (3 — 0.98. For larger /3, we obtained 
smoother curves. 

There are several definitions of the width of 
the fiux tube. Here, we have investigated this in 
terms of the peak radius of the monopole current 
profile, pos- In order to find poff we have interpo- 
lated the on-axis data with a smooth curve. The 
result as a function of r/rg is shown in Fig. ^ 
where we have picked only the data with R = 5 
for all /3 values; pcS seems to grow as increas- 
ing r. However, we make no statement on the 
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Figure 4. The peak radius of the monopole cur- 
rent profile as a function of r/rg. 



behavior whether the width diverges as r — > 00. 
Also, since the rotational invariance is not good 
for small (3 values, we have to check the valid- 
ity of the smooth interpolation carefully. Further 
detailed analyses are in progress 
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